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Abstract. We summarize recent 2D MHD simulations of line-driven 
stellar winds from rotating hot-stars with a dipole magnetic field aligned 
to the star's rotation axis. For moderate to strong fields, much wind 
outflow is initially along closed magnetic loops that nearly corotate as a 
solid body with the underlying star, thus providing a torque that results 
in an effective angular momentum spin-up of the outflowing material. But 
instead of forming the "magnetically torqued disk" (MTD) postulated in 
previous phenemenological analyses, the dynamical simulations here show 
that material trapped near the tops of such closed loops tends either to 
fall back or break out, depending on whether it is below or above the 
Keplerian corotation radius. Overall the results raise serious questions 
about whether magnetic torquing of a wind outflow could naturally result 
in a Keplerian circumstellar disk. However, for very strong fields, it does 
still seem possible to form a centrifugally supported, "magnetically rigid 
disk" (MRD), in which the field not only forces material to maintain 
a rigid-body rotation, but for some extended period also holds it down 
against the outward centrifugal force at the loop tops. We argue that 
such rigid-body disks seem ill-suited to explain the disk emission from 
Be stars, but could provide a quite attractive paradigm for circumstellar 
emission from the magnetically strong Bp and Ap stars. 



1. Background and Context 

Within the broad theme of this conference on "Magnetic Fields in O, B, and 
A Stars" , a key issue is the role of magnetic fields in channeling the line-driven 
stellar wind outflows from such hot, luminous stars. Initial results of full MHD 
simulations for models without stellar rotation (ud-Doula and Owocki 2002) are 
reviewed in these proceedings through the write-up for the talk presented by A. 
ud-Doula. Here we summarize further simulation results for rotating hot-stars 
with a magnetic dipole aligned to the stellar rotation axis. These simulations 
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provide a timely dynamical test of the "Magnetically Torqued Disk" (MTD) 
paradigm, which, as reviewed in the talk by J. Brown in these proceedings, has 
recently been promoted as a central mechanism for producing the disks inferred 
in Be stars (Cassinelli et al. 2002). 

To allow greater emphasis on our MHD results, the broader discussion of 
alternative Be-disk models given in the oral version of the present paper will not 
be reproduced here, as most of this has already been included in the write-up for 
a related talk presented at the recent IAU Symposium 215 on Stellar Rotation 
(Owocki 2003). A key general result of these MHD simulations regards the 
tendency for magnetic torquing of wind material to result in centrifugal mass 
ejection rather a circumstellar disk. 



2. Alfven Radius vs. Keplerian Rotation Radius 

Recent MHD simulations (ud-Doula 2002; ud-Doula and Owocki 2002) indicate 
that the effectiveness of magnetic fields in channeling a stellar wind outflow can 
be characterized by the ratio of the magnetic to wind energy densities 

f]ir) = ^ = r] * (l-R*/ry (1) 

Here 77* = B 2 Rl/(Mv 00 ) defines an overall "magnetic confinement parameter" 
in terms of the strength of the equatorial field B* at the stellar surface radius 
i?*, and the wind terminal momentum Mv^. The latter equality thus isolates 
the radial variation in terms of a magnetic power-law index q (= 3 for a dipole) 
and a velocity index (3 (~ 1 for a standard CAK wind). If, for simplicity, we 
ignore the wind velocity variation (i.e. by taking (3 = 0), we can easily solve for 
an "Alfven radius" t]{Ra) = 1 at which the magnetic and wind energy densities 
are equal 

R A = vl'^R* • (2) 

As shown by simulation results summarized below, this Alfven radius provides 
a reasonable estimate for the maximum radius of closed loops in a wind outflow. 
Moreover, since in rotating models such closed loops tend to keep the outflow 
in rigid-body rotation with the underlying star, it also defines the radius of 
maximum rotational spin-up of the wind azimuthal speed. 

To characterize such rotational effects, let us next define a Keplerian coro- 
tation radius Rk at which rigid-body rotation would yield an equatorial cen- 
trifugal acceleration that just balances the local gravitational acceleration from 
the underlying star, 



GM _ v\ _ V e 2 q R K 
R 2 K Rk Rl 



(3) 



where V eq is the stellar surface rotation speed at the equator. This can be solved 
to yield 

R K = W~ 2/3 R*, (4) 



where W = V eq /V cr it, with V cr u = y/GMjR* the critical rotation speed. 
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Figure 1. The key domains in a parameter plane of magnetic field 
strength (as represented by ^fr\l ~ B*) vs. stellar rotation (as repre- 
sented by the critical rotation fraction W). 



Finally, it is also worth noting here that corotation out to an only slightly 
higher "escape radius" , 

R E = 2 1 I 3 R K = 2 1 / 3 W~ 2/3 R^ , (5) 

would imply an azimuthal speed that equals the local escape speed from the 
star's gravitational field. 

3. MHD Models with Parameters Optimized for Keplerian Spin-up 

The above scalings suggest that a likely necessary condition for propelling out- 
flowing material into a Keplerian disk is to choose a combination of parameters 
for magnetic confinement vs. stellar rotation such that Rk < Ra < Re- In 
the parameter plane of y/rjl ~ B* vs. W defined in figure 1, the required com- 
bination is represented by the relatively narrow gray band. The dark region 
below this represents parameter combinations for which the azimuthal speed at 
the Alfven radius is expected to be sub-Keplerian, while the white region above 
represents cases for which the rotation speed at the Alfven radius should exceed 
the local gravitational escape speed. 

As a sample test case, we focus here on the specific combination r/* = 10 
and W = 1/2, which as shown by the dot in figure 1, lies in the middle of the 
gray domain, and thus should represent an optimal case for magnetic spin-up 
into Keplerian orbit. 

Figures 2 and 3 illustrate results of 2D MHD simulations for this case, using 
the approach and general model assumptions described in ud-Doula and Owocki 
(2002), but now extended to include field-aligned rotation. Figure 2 shows that 
conditions at a time 90 ksec after introduction of the field do superficially re- 
semble a magnetically torqued disk. In the upper panel the darkest regions - 
which represent a near rigid-body rotation with V^jVrigid ~ 1 - include essen- 
tially the entire closed magnetic loop. The log-density grayscale in the lower 
panel shows moreover that magnetic channeling of the outflowing stellar wind 
leads to formation of a dense equatorial compression that might be interpreted 
as the expected circumstellar disk. 



4 



S. Owocki & A. ud-Doula 




Figure 2. Results of a 2D MHD simulation for a model with 77* = 
10 and W = 1/2, shown at at time snapshot 90 ksec after a dipole 
field is introduced into an initially steady-state, unmagnetized, line- 
driven wind. In the lowel panel the grayscale represents the log of the 
mass density, while in the upper panel this represents the ratio of the 
azimuthal speed to the local rigid-body rotation speed, V r ig%d = 
V ro t sin 6 (r/i?*), with white to black corresponing to the range 0.7 to 
1. The curves denote magnetic field lines, and the vertical dashed lines 
indicate the equatorial location of the Keplerian, Alfven, and Escape 
radii defined in eqns. (2), (4), and (5) of the text. 
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Figure 3. Same as lower panel of figure 2, except for time t=390 
ksec, with arrows illustrating the upward and downward flow direction 
of dense material above and below the Keplerian radius. 
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Closer examination shows, however, that most of this equatorial compres- 
sion does not have the appropriate velocity for a stable, stationary, Keplerian 
orbit. Thus in just a few ksec of subsequent evolution, this putative "disk" 
becomes completely disrupted, characterized generally by infall of the material 
in the inner region, i.e. below the Keplerian radius Rk, and by outflow in the 
outer region above this Keplerian radius. Figure 3 illustrates the irregular form 
of the dense compression at an arbitrarily chosen later time (390 ksec from the 
initial start). The arrows emphasize the flow divergence of the dense material 
both downward and upward from the Keplerian radius. This evolution is most 
vividly illustrated through animations, which can be viewed on the web at: 
www.bartol .udel . edu/~owocki/animations/den4wp5etal0 . avi 

We have carried out similar MHD simulations for a moderately extensive set 
of combinations for the rotation and magnetic confinement parameters. In all 
cases we find that any equatorial compressions are dominated by radial inflows 
and/or outflows, with no apparent tendency to form a steady, Keplerian disk. In 
weak magnetic spin- up cases with Ra < Rk, material trapped on closed loops 
generally falls back on the star, much as in the non-rotating models discussed 
in the talk by A. ud-Doula (see also ud-Doula 2002 and ud-Doula and Owocki 
2002). The intermediate magnetic spin-up cases with Rk < Ra < Re (lying 
in the gray area of figure 1) generally show a combination of infall and outflow, 
much as described above. Finally, for the strong magnetic spin-up models with 
Ra > Re the material density tends to build up at the tops the highest loops 
before breaking out in semi-regular epsiodes of discrete "mass ejections". 

Unfortunately technical issues make it difficult to carry out simulations 
with very high magnetic fields and/or very rapid stellar rotation. For rotations 
near the critical speed, the star becomes distorted into an oblate spheroid, thus 
making it difficult to formulate a static atmosphere lower boundary condition 
within our spherically symmetric version of the ZEUS MHD code. For very 
strong magnetic fields (i.e. 77* > 100), the large associated Alfven speed im- 
plies a small Courant time-step within the explicit time-stepping of the Zeus 
MHD code, and so such cases become very computationally expensive to evolve 
through the required few characteristic flow evolution times. Further work is 
thus needed to explore such cases of near-critical rotation and/or very strong 
magnetic confinement. 

4. Magnetically Confined "Rigid-Body" Disks 

In lieu of detailed simulations, it is however possible to infer some likely at- 
tributes of magnetically torqued outflows in the strong field limit, for which 
the magnetic lines-of-force can be idealized as "rigid pipes" that channel the 
wind outflow toward a collision near fixed loop tops. For the case of a ro- 
tation aligned dipole, this effectively concentrates material at the combined 
magnetic/rotational equator. Below the Keplerian radius, such material will 
again tend to fall back to the star, but above this radius, it will be centrifu- 
gally pushed against the loop top, gradually building in mass until it finally 
has a sufficient rotational energy to break open the field. By comparing the 
gravitational-rotational energy of this building disk to the magnetic field en- 
ergy, one can derive an estimate for this "break-out" time. For typical stellar 
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parameters, the ratio of this to a characteristic wind flow time Tfi ow = R*/Voo 
scales roughly as 

Ureak . O.Oll^- 4 / 3 

Tfiou, ~ (r/i?x) 4 - (r/ifo) ' 1 j 

Such a disk would thus have a roughly constant inner edge at the Keplerian 
radius Rk, but a highly variable outer edge. Indeed, having any sort of disk of 
accumulated material that extends much beyond Rk requires that the numerator 
in eqn. (6) exceed unity. In that case, we can estimate a maximum outer cutoff 
radius to be where Tbreak ~ t/Zoum which solves to 

Rout,max ~ ~J=^ • (7) 



A key point here, however, is that such disks would be quite distinct from a 
Keplerian disk, both physically and in terms of likely obervational signatures. In 
particular, they seem likely to be much more variable, characterized by intervals 
with substantial radial flow speeds in the outer disk. Perhaps even more notably, 
their rotation would follow a rigid-body law, ~ r, instead of the Keplerian 
form, v<f, ~ l/xA 7 - 

In terms of a proposed application to Be-star disks, both properties seem 
to run counter to traditional interpretations of Be-disk emission line-profiles 
(Hanuschik 1995), which in many Be stars exhibit features (e.g. "central-quasi 
emissions"; Rivinius et al. 1999) that suggest a quite low upper limit (ca. 10- 
20 km/s) to radial outflow speeds, and moreover seem generally consistent with 
an azimuthal speed that follows a Keplerian law. The problems becomes partic- 
ularly acute for extended rigid-body disks, since these imply more rapid rotation 
speeds than are typically inferred from Be emission-line-widths. For disks with a 
limited radial extent, the differences from Keplerian rotation are less dramatic, 
and as explored in the poster contribution by D. Telfer et al., the resulting 
profiles could indeed be consistent with existing observational analyses. 

However, the long-term V/R variations seen in a substantial fraction of Be 
stars do still seem best explained by a Keplerian disk undergoing a precession 
of elliptical orbits that characterize an one-arm disk oscillation (Savonije and 
Heemskerk 1993; Telting et al. 1994; Savonije 1998). In a rigid-body disk in 
which the individual fluid elements of the disk are tied to the rotation period of 
the star, it is difficult to see how any processes could reproduce the year-long 
timescales of V/R asymmetries that require an associated long-term distinction 
in the physical properties of emitting material in a specific fixed-frame direction. 

On the other hand, such rigid-body magnetic disks may represent a good 
paradigm for the observed circumstellar emission from Bp stars (see, e.g. talk 
write-up by Groote and references therein). The observed ca. 10 kG fields of 
these stars combined with their generally modest mass loss rates M ~ lO~ 9 M /yr 
imply a huge magnetic confinement parameter, of order r/* ~ 10 6 . With the as- 
sociated Alfven radii of Ra ~ 30i?*, even the generally modest rotations of 
stars with, say, W « 0.1, would allow centrifugal disk support above Rk ~ 
10 2 / 3 R* w 4.6-R*, extending perhaps out to R ou t,max ~ 100i?*. Such low-density 
winds are moreover subject to ion runaway instability (Krticka and Kubat 2000, 
2001; Owocki and Puis 2002), and this likely plays an important role in the 
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abundance variations often associated with the circumstellar emission from Bp 
stars. Future efforts should thus explore further such a wind-fed "magnetically 
rigid disk" (MRD) paradigm for Bp stars. 
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